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I present variational solutions of the many-body Schrodinger equation for a two-dimensional 
electron system in strong magnetic fields, which have, at v — 1, the energy about 46% lower than 
the energy of the Laughlin liquid. At v = 2/3, I obtain the energy, about 29% lower than was 
reported so far. 

PACS numbers: 73.40.H 

The need to explain the integer |l] and the fractional g quantum Hall effects || (QHE) resulted in tremendous the- 
oretical efforts |J-|2l[] to understand the nature of the ground state of a system of strongly interacting two-dimensional 
(2D) electrons in strong magnetic fields. The modern, widely accepted explanation of the fractional QHE is based 
on the theory of Laughlin ||, who showed that a special quantum state - an incompressible quantum liquid with 
fractionally charged quasiparticlcs |22-27| - has the lowest ground state energy at the Landau- level filling factors 



v = 1/m, m — 1,3,5. In this Letter I present an infinite set of variational many-body wave functions for a 2D 
electron system (ES) in strong magnetic fields B, and an exact analytical method of calculating the energy of the 
new states. The new wave functions describe the system at all Landau- level filling factors v < 1, continuously as a 
function of magnetic field. In certain intervals of v (approximately, at v > 1/2) the calculated upper limits to the 
ground-state energy are substantially lower than was reported up till now P Jl^JT^] . 

The main theoretical efforts so far have been concentrated on the region of the very strong (y < 1) magnetic fields 
[p]-|l9||. The nature of the ground state of a completely filled lowest Landau level (y = 1) seemed to be well understood: 
the Hartree-Fock many-body wave function 

* if F = ^det|^(r l )|, Lj = 0,1, ... ,N — 1, (1) 

- a Slater determinant built from the single-particle lowest-Landau-level states 

^(r) = -^4=exp(-zz72), (2) 

- was proposed about 20 years ago JSTJl and was considered to be the cornerstone of the theory [here z — (x + iy)/\ is 
a complex coordinate of an electron, A 2 = 2hc/eB, L = 0, 1, 2, ... is the angular momentum quantum number, and h, 
c, and e are the Planck constant, velocity of light and the electron charge, respectively]. The state (|l]) is a special case 
of the Laughlin liquid. It is characterized by the uniform electron density n s = 1/ttX 2 and the energy per particle || 

e H F = -tt/2 = -1.5708 (3) 

in the ^-independent energy units e 2 yfnl. 

The wave function (|lf) describes electrons rotating around the same center (arbitrarily placed at the origin) with 
different angular momenta. However, in an infinite and uniform 2DES all points of the 2D plane are equivalent and 
no one of them can be considered as a special point, which electrons would prefer to rotate around. I consider, as 
trial functions for the ground state of the 2DES, an infinite set of many-body wave functions 

* L = -^=det|4|, (4) 
ViV! 



4 = XLiTi,*,) = - R i ) e -^-(BxR,)/^ ) (5) 

which describe a system of electrons rotating around different points Rj with the same angular momentum. Here 
<f>Q = hc/e is the flux quantum, and L = 0,1,2,.... The points Rj are chosen in view of minimizing the Coulomb 
energy of the system, and in a pure 2DES (without disorder) coincide with points of a triangular (or a square) lattice, 
uniformly distributed over the 2D plane with the average density n s (n s a 2 — 2/V3 and 1 for the triangular and the 
square lattice, respectively, a is the lattice constant) . All the states VP l are the eigen functions of the kinetic energy 
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operator, with the energy hu> c /2 and the angular momentum L per particle. The state ^l=o is the Wigner crystal 
state 0. The states ^l>o have been discussed in the literature ^8|, but only in Hartree approximation, when the 
neighbour single-particle states tpL do not overlap. As will be seen below, the most interesting physics takes place 
just in the opposite case of a strong overlap of the states ipL centered at different lattice points Rj. 

In order to calculate expectation values of different physical parameters with the wave functions (|4|), I derive exact 
analytical formulas for the norm 

(* L \* L )=det\S\ ) (6) 

and matrix elements of arbitrary single-particle, l\Hi\^ l) , and two-particle, l\H2\*& l) , operators. The norm 
(|J) and the matrix element of Coulomb energy, 

(*l\V c \* l ) = \ ]T(-iy Wsgn(* - fc) S gn(j - l)V^ det \S\ r ™*^ v ( 7 ) 

ijkl 

are expressed via the matrix S with the elements 

S ij = S% L = { X L(r,-R i )\XL(r,R j )), (8) 

and the Coulomb matrix elements = (XaiXbk\ e2 \ r a ~ r b\~ 1 \XajXu) ■ The sums in (^) are taken over N 3> 1 lattice 
points. Both the overlap integrals S^ L and the matrix elements V^j can be analytically calculated, so that I get, 
finally, exact closed-form analytical expressions for the energy el of the system (per particle) in the states Vf^. The 
final result can be presented in the form = + ef, c - The Hartree contribution eff contains the sum J2i<j ^ujj an d 
the energy of the uniform positive background. The exchange-correlation energy contains two-site exchange terms 
Vijji, as well as three-site (Vuki, Vijki) and four-site (Vyjy) correlations. 

The number N of lattice points involved in calculation of sums @) should, ideally, be infinite. I calculate the Hartree 
energy exactly, in the thermodynamic limit N = oo, and the exchange-correlation energy e : f j c (N) approximately, at 
a finite number of lattice points. It turns out that the contribution e| / c (iV) is negative and grows with N in its absolute 
value, Figures ||a and[2| The calculated energy el(N) is thus an upper limit to the true energy cl, £l < €l(N). 

Figure [j] exhibits B-dependencies of the calculated energies el (N) at a number of different L and N in a system 
with a triangular lattice of points Rj. As expected, in strong magnetic fields the states with larger L have larger 
energy, due to a larger Hartree contribution. When B decreases, however, the single-particle wave functions rotating 
around different lattice points begin to overlap, and exchange-correlation interaction significantly reduces the energy. 
Due to the overlap of electron wave functions and electron interference effects, the energy of the system in the states 
has a complex oscillating dependence on magnetic field. As a consequence, in different intervals of B the ground 
state of the system is realized by states with different L. The gap between the ground and the first excited state also 
varies with B and vanishes in points, where two different L-states have the same energy. It seems likely that, maxima 
in the diagonal resistivity of the system are observed in these, zero-gap, points, while the Hall-resistivity plateaus are 
seen in the regions of large gaps, when the system is in the same quantum state. This view on the integer and the 
fractional QHEs merits detailed consideration and requires further development of the theory. 

At v = 1, the state ^>l=3 has the lowest energy. For this state I calculate e^ =3 {N) with inclusion of up to N = 187 
lattice points, Figures [l] and g. Thus calculated upper limit to the total energy per particle, 

e L=3 {u - 1) < e L= ^ N=l67 {u = 1) = -2.302, (9) 

is more than 46% lower than the energy (||) of the Hartree-Fock state (|l|) . Note that the energy (||) is even lower than 
the energy of the classical Wigner crystal |2{| . 

The state ^l=3 at v = 1 is characterized by a very strong overlap of the single-particle wave functions of neighbour 
electrons, Figure || As a consequence, the density of electrons n e (r) in this state is uniform (left panel). The wave 
function \& l=z{v = 1) thus describes a quantum liquid with strong exchange-correlation interaction of a very large 
number of electrons. [For comparison, the right panel shows the electron density in the same state {L = 3) at v = I /5 
(the state v Pl=3 is one of the excited states at v = 1/5, Figure [l]b). Qualitative difference between the liquid-type 
state at v = 1 and a solid- type state at v = 1/5 is obvious]. 

The wave functions \1/ l have a very large variational freedom. Apart from a possibility to vary the angular momen- 
tum L, one can consider configurations with a different (e.g. square) symmetry of lattice points Rj. For instance, I 
have observed that, at v = 2/3 the square lattice configuration gives the energy, lower than the triangular one. The 
upper limit to the total energy per particle at v — 2/3, calculated with L = 5 and N — 149 in the square lattice 
configuration (the filled diamond in Figure [l]b), 
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e L=b {v = 2/3) < e L= ^ N= Mu = 2/3) = -2.050, 



(10) 



is about 29% lower than was reported so far |l7j and also lower than the energy of the classical Wigner crystal. 

Can the energy of the states $i be comparable with or lower than that of the Laughlin liquid at v = 1/3 and 
v = 1/5? This important question remains open. Calculations performed so far for L < 3 and N < 91 gave the 
energy higher than that of the Laughlin liquid at these values of v. However, Figure [l] demonstrates two important 
tendencies: i) the downward cusps in dependencies cl(B) are being shifted to the right with growing L; and ii) for 
given L and v the energy £l(./V) substantially decreases with the growth of N. Calculations with other L and larger 
N may lead to a reduction of energy ex below the energy of the Wigner-crystal state el=o also at v < 1/3. 

The new approach presented here enables one to analytically calculate any physical value characterizing the system 
at all v , continuously as a junction of magnetic field. It naturally describes a transition to the Wigner crystal state 
at v — > 0, and can be straightforwardly generalized to the case of non-spin-polarized systems and hence to the region 
of higher Landau levels v > 1. It also makes possible to study the interplay between disorder and electron-electron 
correlations in the integer and the fractional QHE (by varying configurations of points Rj which are influenced by 
disorder in a real system |30| ) . Due to a large variational freedom of wave functions it opens up wide possibilities 
to search for better approximations to the ground state at all values of v. 

This work was supported by the Max Planck Society and the Graduiertenkolleg Komplexitat in Festkorpern, 
University of Rcgcnsburg, Germany. I also thank Peter Fulde and Ulrich Rossler for support of this work, Nadejda 
Savostianova for numerous helpful discussions, and Vladimir Volkov for useful comments. 
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FIG. 1. Calculated energies 6l(N) — +ef, c (iV) per particle (in units e 2 y / n7) as a function of magnetic field (of the inverse 
filling factor v^ 1 = B/n s cj>o): a, at L = 3 and different numbers of lattice points N; b, in the states $ l with a few lowest L. 
The filled diamond in b corresponds to a square lattice of points Rj , all other points and curves - to a triangular lattice. The 
Hartree contribution e L (the uppermost curve in a, thin curves in b) is calculated exactly in the thermodynamic limit N = co. 
For comparison, the energy of the classical Wigner crystal and of the Laughlin liquid |9|Jf^] are also shown. 
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FIG. 2. Absolute value of the exchange-correlation energy per particle, —ej^{N), in units e 2 ^/n^, as a function of the number 
of (triangular) lattice points N, for the state ^l=3 at v = 1. 




FIG. 3. Calculated normalized density of electrons n e (r)/n s in the state ^l=3 at v = 1 (left panel) and v = 1/5 (right 
panel). Up to iV = 301 lattice points have been included in this calculation. Insets shows an overlap of electron rings centered 
at sites of the triangular lattice Rj . 
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